from scipy import stats
import numpy as np
import matplotlib.pyplot as plt
n=6
p=0.30
X=np.arange(0, n+1,1)
pList = stats.binom.pmf(X, n, p)
plt.plot(X, pList, marker='o',linestyle='None')
plt.rcParams['font.sans-serif']=['SimHei']
plt.vlines(X, 0, pList)
plt.xlabel('随机变量：购买服装的顾客人数')#x轴文本
plt.ylabel('概率')#y轴文本
plt.title('二项分布：n=%i,p=%.2f' % (n,p))#标题
plt.show()#显示图形
print("恰有四个顾客购买的概率为：%f"%(stats.binom.pmf(4, n, p)))
print("购买的顾客不超过一半的的概率为：%f"%(1-(stats.binom.pmf(4, n, p)+stats.binom.pmf(5, n, p)+stats.binom.pmf(6, n, p))))
print("至少有一个顾客购买的概率为：%f"%(1-stats.binom.pmf(0, n, p)))